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Part of the CMB polarization signal in the direction of galaxy clusters is produced by Thomson 
scattering of the CMB temperature quadrupole. In principle this allows measurement of the CMB 
power spectrum harmonic C2{z) with higher accuracy (at 2 > 0) than the cosmic variance limit 
imposed by sample variance on one CMB sky. However the observed signals are statistically cor- 
related if the comoving separation between the clusters is small enough. Thus one cannot reduce 
the sample variance by more than roughly the number of separate regions available which produce 
uncorrelated signals, as first pointed out by Kamionkowski and Loeb. In this paper we analyze in 
detail the procedure outlined by Kamionkowski and Loeb, computing the correlation of the polar- 
ization signals by considering the variation of the spherical harmonic expansion coefficients of the 
temperature anisotropy on our past light cone. 

Given a hypothetical set of Stokes parameter measurements of the CMB polarization in the 
directions of galaxy clusters, distributed at random on a given redshift shell, we show how to 
^1^' construct an estimator of the angular power spectrum harmonic C2 at that redshift. We then 

CD ■ compare the variance of this estimator with the cosmic variance of the CMB multipole on our sky 

which probes the same scale. We find that in fact the cosmic variance is not reducible below the 
QQ ' single sky CMB value using the cluster method. Thus this method is not likely to be of use for 

reconstruction of the primordial power spectrum. However the method does yield a measurement 
, of C2 as a function of redshift with increasing accuracy at higher redshift, and thus potentially a 

^ ' probe of the mechanism which may have suppressed the quadrupole. 

fT^ ' We also examine to what extent the redshift dependence of C2 can be used to probe the time 

, changing potential anisotropy as the universe evolves into the vacuum dominated phase (the late- 

T-H ■ time integrated Sachs- Wolfe effect). We find that this effect is not observable in the time dependence 

' of C2 since it is swamped by cosmic variance, but there is an observable signature in the correlation 

, functions of the Stokes parameters. 

^ . 

O ■ PACS numbers: 98.80.Es,95.30.Gv,98.70.Vc 

o ■ 

^ ; I. INTRODUCTION 

C/3 . 

, The CMB radiation incident on galaxy clusters has an intrinsic intensity quadrupole Q2 created by inhomogeneity 
' at the surface of last scattering. Thomson scattering of the CMB in a galaxy cluster with typical line of sight optical 
k>( I depth rc generates polarization of order Q2Tc- Thus a measurement of this polarization signal would allow an estimate 
' of the CMB quadrupole at non-zero redshift. This is of interest because it would potentially allow us to get around the 
restriction of cosmic variance. To elaborate, at z = we only have one CMB sky to observe, with {21 + 1) independent 
real data points for each spherical harmonic mode of the CMB on our sky, to compare with the ensemble average 
prediction of the variance. There is thus an intrinsic fractional sample variance of the harmonic Ci of 2/ {21 + 1) (see 
section iflTI) . which severely limits comparison with the ensemble averaged theory at low I. This restriction limits the 
accuracy of measurements of the primordial power spectrum on the largest scales. The theoretical predictions thus 
obtained for the CMB power spectra are fundamentally limited by this sample variance, commonly termed the cosmic 
variance. 

Thomson scattering of the I ~ 2 part of the CMB anisotropy in a cluster generates a secondary polarization 
anisotropy which depends on the spherical harmonic components a2m as seen by a (hypothetical) observer at the 
cluster. Since this polarization signal produced by a cluster is sensitive to the density perturbations on a last scat- 
tering surface different to our own, this in principle allows one to make more accurate comparison to the theoretical 
predictions for CMB angular power spectra at low / than allowed by the cosmic variance limit. However the observed 
signals are correlated if the comoving separation between the clusters is small, and many strongly correlated signals 
arc no more useful for reducing the sample variance than one signal. The variance in the estimated quadrupole can 
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be reduced by roughly the number of regions available which produce uncorrelated signals. 

This method of using the CMB polarization signal produced by galaxy clusters to get around cosmic variance limits 
was first pointed out by Kamionkowski and Loeb (we usually refer to it as the "cluster method" in what follows). 
However Kamionkowski and Loeb did not actually compute the correlation of the cluster signals in a particular 
cosmological model in their paper, and did not therefore demonstrate explicitly that the cosmic variance is reduced 
with a given set of clusters, nor did they develop any formalism for converting measurements of the Stokes parameters 
of the CMB to statistical estimators which get around cosmic variance. In [J, estimators of C2(t) were constructed 
(taking into account the kinematic SZ contamination of the polarization signal also), but the effect of statistical 
variation in the polarization signal on the estimator variance was not included. This variation was considered by |3l , 
but they computed the variation of the quadrupole as an expansion in small cluster separations, and their analysis 
is not applicable to a general set of clusters in arbitrary locations. In this paper we compute the correlation of 
the cluster signals in the case of an idealized set of measurements from clusters distributed in random directions on 
a given redshift shell. We describe an explicit procedure for carrying out the program outlined in and study 
how the correlations die off as clusters of increasingly high redshift are used. Information about the correlation 
of the polarization signals is contained in the generalized correlation functions of the CMB temperature anisotropy 
coefficients, (a;m(a;, r)a;*^,(a;', r')), which contain all of the statistical information (assuming Gaussianity) about the 
variation of the aim coefficients as the observation point and associated last scattering surface change. With these 
functions, we can derive an estimator for C2{z) in terms of Stokes parameters, and find its variance. 

We should clarify exactly what we mean by reducing cosmic variance. It is true that the polarization signals provide 
an estimator C2(t) of the remote quadrupole at a given redshift which has a smaller fractional cosmic variance than 
the local quadrupole. However, this is not the most useful comparison, since this estimator probes smaller physical 
scales than the local quadrupole. Cosmologists already have estimators of the power on these scales, namely the 
WMAP angular power spectrum harmonics C; with I > 2. So the interesting question to ask is whether (r) has 
smaller cosmic variance than the CMB multipole on our sky that probes the same physical scale. This determines 
whether or not the cluster technique is capable in principle of providing a better reconstruction of the primordial 
potential than the WMAP data. The results are presented in jvl 

We now outline the organization of the paper. In ^llll we derive the two-point generalized correlation functions of 
the spherical harmonic coefficients, assuming a Gaussian primordial perturbation spectrum. In illlll we discuss the 
CMB transfer functions used to compute the generalized correlation functions, and examine the time dependence of 
C2(t). In illVI we derive expressions for the Stokes parameters Q, U (defined in an appropriate all-sky basis) of the 
CMB radiation scattered into the line-of-sight by the cluster gas, in terms of the local aim at the cluster, for a general 
line-of-sight. Note that in this section we found it convenient to use the "density matrix" formalism for polarization 
calculations 0, 0| , which is outlined in Appendix ^ Then in ^3 we consider the the statistical variation of these 
Stokes parameters with the comoving position of the cluster. We construct a simple estimator for C2(t) and compute 
its variance for a number of simulated sets of clusters. In WII there is a discussion and summary. Note that we restrict 
the discussion to the case of a flat FRW universe throughout, for simplicity. 



II. GENERALIZED CORRELATION FUNCTIONS 



The observed CMB sky and associated power spectrum changes with the comoving spatial position and redshift 
of the observer. The statistical variation with position may be characterized completely with a set of correlation 
functions of the temperature perturbation. 

The fractional CMB temperature perturbation A(a;, n, r) AT/Tcmb is a function of the comoving spatial position 
of the observer x, the conformal time of observation r, and the direction of the line of sight n. In order to expand 
the CMB temperature field in spherical harmonics aim{x,T) we must define a polar coordinate system everywhere. 
In flat space, it is simplest to use the convention that the polar axis is taken to be the same at each point (in a curved 
space, this would not be possible and more care would be needed). 

When we need to specify coordinates for n we use spherical polars {9, 0), with the polar axis aligned with the 
direction and the = plane normal to the By direction. The fractional temperature perturbation is expanded as 

oo / 

A(a;, n,T)=^ ^ airn{x,T)Yim{n) ■ (1) 

;=o m=-l 

The two point correlation function of the temperature anisotropy is 

(A(a;,n,r)A(a;^n^TO) =^^G„(a^,T,a;',T')Fz„(n)y;„,(n') , (2) 

Im I'm' 
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FIG. 1: Comoving coordinate system for generalized CMB correlation function, in flat space. All points in the plane shown 
lie on the past light cone of the observer at O. The outermost circle indicates our last scattering surface. The circles centered 
on the observers at positions x and x' (and conformal times r = ro — \x\ and t' = to — \x\ respectively) indicate their last 
scattering surfaces, which are smaller since recombination occurred in the less distant past according to them. Note that any 
orientation of the points x and x' in space can be rotated into this plane. 



where we have defined the generahzed CMB correlation functions 

Cimi'm'{x,T,x' ,t') EE (fl/™ (a;, r)a;* (aj' , r')) . (3) 

The set of functions CimVm'{x,T,x' ^t') form the covariance function of the Gaussian random process from which 
the photon distribution function, defined at all points in space and observed in all directions, is sampled. 

Given a set of cosmological parameters, the generalized correlation functions may be computed using the photon 
transfer function which describes the physics of the propagation of CMB photons from the last scattering surface to 
the cluster. The generalized correlation functions manifestly obey some simple symmetry relations: 

Cimi'm'{x' ,x) = C;*,„,;„j(a;, a;') . (4) 

We will be interested only in the case where both sets of spacetimc coordinates lie on our past light cone, in order 
that we arc computing only quantities which arc directly observable. Choosing ourselves to be at the spatial origin 
of the comoving coordinate system, at conformal time tq (the age of the universe in conformal time), the events at 
X, x' occurred at conformal times r = tq — |a;|, t' = tq — |a;'| respectively. The correlation function may therefore be 
written as a function of spatial variables only, C;mi'm'(^: ^')- The coordinate system is illustrated in Figure ^ Note 
that the last scattering surfaces of observers at z > are spheres which are tangential to the last scattering sphere of 
an observer at z = 0. 

Now we briefly discuss the statistical properties of the coefficients aim and review the usual definition of cosmic 
variance. The aim coefficients at any given point are all independent, but there are spatial and temporal correlations 
between any pair of coefficients at different points. Assuming Gaussianity of the primordial perturbations, it follows 
that all of the various A'^-point joint probability distribution functions for the aim at separate points are (complex) 
multivariate Gaussians, with covariance matrix R given by (note that each index labels both the set of values Zi,mi 
and also the point Xi in three dimensional space): 

Rij = {0'li7ni{Xi)ai.m-{Xj)) = Cl-miljmjiXi, Xj) . (5) 

Thus given a cosmological model wc have the p.d.f of the ensemble from which the aimix) arc drawn, and the 
associated ensemble average angular power spectrum harmonics Ci(t). Then given a set of observations o°^(a;), these 
ensemble average predictions are typically compared to the observed quantities C°(a;) = \a°^{x)\'^ / (21 + 1) 

(clearly (Cf (a;)) = Q(r)). 
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On the sky of an observer at time r, the mean square difference between the observed CMB angular power spectrum 
C°{x) and the ensemble average theoretical power spectrum Ci{t) is characterized by the cosmic variance 

([C^ix) - Q{r)r) = {Ctixf) QirY . (6) 

Expanding we obtain 

(Cnx?) = ^(2/ + l)-2(|aL(a:)naL'(^)r) • (7) 

mm' 

To evaluate the right hand side, we need the ensemble average of the product of four a;m's. This follows from 
Gaussianity: 



(a;,„j(a;,T)a,*„^(a;,r)ai'„3(a?',r')a^,,„^(a;',r')) = (5,„im2'5m3mjCi(T) Ci>{t') 

+ Chnil'ynA^^^ )^lm.2l' m-ii.^ i ^ ) + C'imii',-m3 (^i ^ )Q*n2i' .-m.j (^J ^ ) • (8) 



Now setting I' ^ I, x = x', and using Cimi'm'{x,x) = 5ui5mm'Ci{T)^ we obtain the familiar expression for the cosmic 
variance associated with each harmonic, 

{Ct{xf)-C,{Tf = ^Qirf . (9) 

This quantity captures how much we can expect the measured power spectra to differ from the ensemble average. 
Note that this expression also follows from the fact that C° is the sum of squares of independent identical Gaussian 
random variables, and is therefore distributed as a scaled X21+1 random variable. 

Related results concerning the spatial correlations of various quantities can be derived similarly. The ensemble av- 
erage difference between the aim measured by separated observers on our past light cone, for example, is characterized 
by 

{\ai,m'ix')-airnix)\'') = Ci^ {t' ) + Ci{t) ~ 2Rc Qml' m' {x , x') , (10) 

and the ensemble average difference between the angular power spectrum harmonics as a function of spatial separation 
is characterized by 



{[cnx)-cr,ix')r) = {crixf) + {cr,{xr)-2{cr{x)c?,ix')) 



- t2l+mi' + l) ^ [\Clml'm'{x,x')\' + \Cl,nl',-m'{x,x')\'] . (11) 
mm' 

We now derive the relationship between the generalized correlation functions and the CMB transfer function. In 
a flat FRW space, the temperature anisotropy may be Fourier expanded in comoving wavcnumber k on the three- 
dimensional hyper-surface of constant r, E,-, 

A{x, n,T)= J <fk e^^ ' =^A(fc, n; E^) , (12) 

where A(fe,n; E^.) is the Fourier transform associated with this hyper-surface only 

Since each Fourier mode corresponds in the case of scalar perturbations to a plane wave perturbation which has 
azimuthal symmetry about fc, A(fc, n; E^) depends only on fe • n and |fc| and may therefore be expanded in Legendre 
polynomials: 

oo 

A(fc,n;E,)=^(-i)'(2/ + l)Az(fc,T)Pi(fc-") • (13) 

The {—if is included by convention to be consistent with most other authors. It is convenient to use the addition 
theorem at this point to express the Legendre polynomials in spherical harmonics, giving 

/oo I 
e^^ • ^^(-i)' AKfc,T) ^ Y:„Xk)Yim{n) ■ (14) 

1=0 m.=~l 
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Employing the orthogonality relation for spherical harmonics J dft ^^(w)!^/™/ (n) = Si'iSm'm (where dfl is the integral 
over solid angle elements centered about direction n) yields 



(15) 

The correlation function may now be written as 

Ci„a'.n'{x.x') = (-z)'-'' {Anf j d^kd'k'e^^-^e-^^' 

{Ai{k,T)A*{k\T'))Y*Jk)Y,rn'ik) . (16) 

Now since the Boltzmann equation which governs the evolution of the CMB anisotropy A((fc, r) is independent of 
k (in linear theory), the k dependence comes entirely from the initial conditions, and we may write A;(fc, r) in terms 
of the CMB transfer function Ai{k, r) which is defined by 0: 

Ai{k,T)^Mk)Mk,r) , (17) 

where 4>i{k) is the initial potential perturbation and Ai(fc,r) is real By the assumption of translational invariance 
4>i{k) has a two-point correlation function of the form 

{Ukm{k'))^P^{k)5^{k~k') , (18) 

where P(j,{k) is the power spectrum of the primordial (post-inflationary) gravitational potential fluctuations. Then we 
may write 

(A,(fc,T)A;,(fc',T')) ^ Aiik,T)A,ik',T')P4k)S^k-k') . (19) 

If we evaluate Cimi'm' {x, x') for x ~ x' (and r = r'), the covariance matrix is diagonal and the familiar orthogonality 
relation follows 

Ciml'ra'{x,x) = {'iTrfSvAa-mj k^dk Af{k,T)P^{k) = (5;'i<5rn'mC; (t) , (20) 

where Ci{t) = {\aim{x,T)\'^) is the ensemble average of the I harmonic of the CMB power spectrum according to an 
observer at conformal time r. Thus at any given point all of the aim are independent random variables. Using the 
addition theorem, we obtain the usual real-space angular correlation function, at any epoch 

{A{x,n,T)A{x,h\T)) = j-J2i^l + ^)Ciir)Piin-n') ■ (21) 
With a; 7^ a;', T 7^ r' we obtain a more general expression for the correlation function: 

Ai{k,T)Ai,{k,T') P^{k) Yi*^{k) Yi,m'{k) . (22) 

The symmetries stated in equation Q may be verifled from this expression. 

We may perform the angular part of the d^k integral in Eqn. H22|) by expanding the plane wave piece in spherical 
Bessel functions (valid in flat space). In a flat FRW space, we are free to manipulate the comoving 3- vectors of events 
as if we were dealing with position vectors in Euclidean space (see e.g. 0, P-71, Eqn. (3.19)). Thus we may deflne 
direction vectors Aa; = a; — a;', Ax = (a; — a;')/|a; — a;'|, and expand 

oo 

exp(ifc-Aa;) = ^{21" + Pin{k ■ Ax)jin{k\Ax\) 

l"=0 

oo /" 

= 4nJ2rj,,{k\Ax\) yi'',n"{k)Yi„,n"{^&) , (23) 

l"=0 m"=-l" 
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where we used the addition theorem to separate the k and A;r dependence. Then the correlation function becomes 
Ci^v^,{x,x') = {Anf J edk Ai{k,T)Ai,{k,T') P4,{k) H)'-'' (-1)™ 

l"m" •' 

The angular integral of the product of three spherical harmonics is expressible in terms of the Wigner 3j symbols (see 
e.g. i), 



V (J.\ V (i.\v* I i2l + mi' + l)i2l" + l) ( I V V'\( I V I" \ 
d^u YUk) Y,m'{k)Y,,^.{k) = y [o )[m m' m" ) ' ^^5) 

The 3j symbols are non-zero only if m + m' ~ m" , and I, I' , I" satisfy the triangle condition that /" be equal to one 
of / + Z', ^ + /' — 1, • • • , |/ — The sum over I" therefore reduces to a finite sum. We have finally 
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' I V l"\f I V I" 



X 

i"=\i~i'\ 



E ( )( jvTO ^'^".H|A.|,r,.')^.^.^-™(A*) , (26) 



where all of the physical information is contained in the kernel 

AV,i"(|Aa;|,r,r') = J k^dk Ai{k,T)Ai,{k,T') P^{k) ji„{k\Ax\) . (27) 

III. TRANSFER FUNCTIONS 

To compute CimVm' for a given cosmological model we need the CMB transfer function A/(A:,r). On the large 
angular scales accessible via the polarization technique, the only significant effects responsible for the temperature 
anisotropy which need to be included in the transfer function are the Sachs- Wolfe (SW) and integrated Sachs- Wolfe 
(ISW) effects 0. The SW effect is the anisotropy due to the gravitational potential fluctuations on the last scattering 
surface, and the associated time dilation effect. The ISW effect arises because, at late times, as the universe is 
making the transition from the matter dominated phase into the vacuum dominated phase, the fluctuations in the 
gravitational potential - on scales still in the linear regime - are still evolving with redshift. As photons fall into and 
climb out of this time changing potential they arc red-shifted and thus a temperature anisotropy is generated. 

We first consider the transfer function of the SW effect. This is computed by ignoring the physics on scales 
comparable to the acoustic horizon at the time of recombination, and retaining only the large scale effects. In this 
limit, the anisotropy is produced solely by the variation in potential <j) (and the consequent gravitational redshift and 
time dilation effects on the photons) and photon density 5^ across the last scattering shell, ignoring the small scale 
acoustic waves which give rise to the acoustic peaks in the angular power spectrum. Using the line-of-sight integration 
method the SW temperature anisotropy is given in real space by 



A{x,n,T) = / dx'C(T-x') 
Jo 



i(5^(x'ri) + 0(x'n) 



(28) 



Here x' is the comoving distance measured along the past light cone of the observer at (a;, r), in the direction n. The 
visibility function is defined by Qir) = e""^*^'"^', with Thomson optical depth tc(t) = dx'a{T — x')''^e{T — xO'^t 
(here and elsewhere a dot means a derivative with respect to conformal time r). 

In the Sachs- Wolfe approximation (valid on scales much larger than the acoustic horizon), and assuming adiabatic 
initial conditions, a perturbative analysis of the equations of motion shows that 5^ = — \(t>, and that in Fourier space 
the evolution of the potential is given by (j){k,T) — j^(f)i{k) (see, for example, 0,0|). The factor of 9/10 accounts 
for the evolution of the transfer function between radiation and matter domination (in the case of adiabatic initial 
conditions). Decomposing the plane waves (working here in flat space) into spherical waves, we obtain 

Mk,r) = I. J\x' ({r ~ X') Jiikx') ■ (29) 
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The visibility function contains the physics of recombination. It rises rapidly during recombination from to 1, 
with derivative sharply peaked about the time of recombination, The effect of the finite thickness of the last 

scattering shell can only influence the radiation field on rather small scales, so for times well after recombination, and 
for low I, we may assume that recombination occurred instantaneously at time r,.. In this limit we may take C to be 



a delta function centered on r^, and the transfer function reduces to 



(30) 



Taking the usual power law form P,j,{k) — Ak"^ ^ (n = 1 gives the scale-invariant Harrison-Zeldovich spectrum) with 
an arbitrary amplitude A (with dimensions of 

{c/ Hq)"' ), with the transfer function in equation (|30|l the integral in equation H2U|) may be done analytically, yielding 



the well known Sachs- Wolfe expression for the angular power spectrum at low / and n < 3 (see for example |3) 



Ciir) 



A(4^)2 J^i-^ k^dk k 
.„-4r(3 



n-4-2 



A{4:TTy 



TTZ 



{r-TrY 



(31) 



Note that if n = 1, this expression has no time dependence. This is a manifestation of the scale invariance property 
of the n = 1 case. 

In the general case including the ISW effect (and assuming a flat universe), the CMB transfer function is given by 
a generalization of Eqn. ()28|l . the line of sight integral: 



A(x,n,T) = 



dx C(t - x) 



S^{x'n) + (l){x'n) 



dx' C(r - x') mx'n) 



(32) 



In linear theory, the growth of the amplitude of the potential perturbations is governed by the growth function D-^- (t) 
of the dark matter perturbation The evolution of the potential perturbation in the adiabatic case is then given by 
(p{k,T) = jQ(t>i{k)D^{T) /a{T). In the case of a flat universe with only non-relativistic matter and vacuum energy, the 
solution for the growth function, normalized to = a at early times, has the simple form |l4j : 

ra{T} 



X'^'ia') da' 



where 



Xia) 



In the instantaneous recombination approximation, this leads to the following CMB transfer function 



A,(/c,r) = —ji [fc(r-T,)] + 



dx' ji{kx') 



d_D^ 
dr a 



(33) 



(34) 



-x' 



x!- 



where the time derivative in the integrand is evaluated at time t 

In the kernel Kunn, the product of two transfer functions appears. Thus there are three terms, a contribution 
entirely from the SW effect, an "interference" contribution from both the SW and ISW effects, and a contribution 
entirely from the ISW effect. Since the ISW part of the transfer function is usually negative, the interference term 
tends to cancel the third term. This is illustrated in Fig. |21 which shows the redshift dependence of the CMB harmonic 
6*2 in a ACDM model with (ilA,r2m) = (0.7,0.3), and various values of the power spectral index n. 



IV. SCATTERING OF CMB QUADRUPOLE 



In this section we consider the generation of polarization by scattering of the intrinsic CMB quadrupolc from 
electrons in a galaxy cluster, which we idealize as a concentrated clump of stationary electrons with a Thomson 
optical depth tc, at a specific angular point on the sky. 
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FIG. 2: A plot of C2{z), including both the SW and ISW effects, normalized to the value of the SW contribution at z = 0. 
The growth function was computed in a ACDM model with (Oa, ^m) = (0.7, 0.3). The upper solid black curve is for the case 
of power spectral index n = 1, and the lower solid black curve for n — 0.9. The dashed upper and lower curves curves show 
the contributions to the n = 1 case from the ISW-ISW term and the SW-ISW interference term respectively. These two tend 
to cancel. (Note that the interference term is negative, and its magnitude is plotted here). 

The Stokes parameters of the radiation scattered into the line of sight to the cluster are functions of the quadrupole 
anisotropy in the local CMB radiation field at the cluster. This is characterized by the coefficients a2m{x) of the 
spherical harmonic expansion of the fractional temperature anisotropy of the radiation field, which are functions of 
the spatial position of the cluster in comoving coordinates denoted x (see for a description of this coordinate 
system). The direction vector of the line of sight from the observer at z = to the cluster is x. The coordinate 
system used is illustrated in Fig. 13 

In terms of the general set of coefficients aim{x), we may write the brightness temperature of the incident CMB 
radiation field at the cluster as a function of the direction vector n of the incoming photon as viewed from the cluster: 

/(n,a;) -TcMB(r)^ {x)Yijn(n) , (35) 

Irn 

where t = tq — \x\ is the conformal time of the scattering events. Since the primary anisotropy has a blackbody 
spectrum, there is no frequency dependence in /(n, x). 

We assume that the incident CMB radiation is unpolarized, which is sufficient to compute the lowest order polar- 
ization signal generated by the quadrupole anisotropy (there are relativistic corrections to the effect described here 
in the case of a cluster with a pecufiar velocity with respect to the CMB, as discussed by which turn out to be 
negligible) . 

The brightness temperature polarization matrix ~ (EiEj) of the radiation scattered into the line of sight is 
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FIG. 3: Illustrates the coordinate system used to describe the generation of polarization by Thomson scattering. The CMB 
incident on a cluster at comoving position x, which we approximate as unpolarized, is Thomson scattered and re-radiated by 
free electrons in the cluster, producing partially polarized radiation scattered into the line of sight. At the observer position 
(redshift 2 = 0) the radiation is decomposed into Stokes parameters with the polarization basis vectors ei, 62 indicated (defined 
in equation 1471 . The CMB radiation incident on the cluster at x is decomposed into spherical harmonics defined with respect 
to polar coordinates 6' and 0'. Only the 1 = 2 harmonics of the incident radiation field generate polarization. 



given in the Thomson limit by the following equation: 

hjix) = {5ik - XiXk) {6ji - XjXi) J dQ' /(n, x) {Ski - fikni) , (36) 

where dft' is the solid angle element about the n direction. A self-contained derivation of this form of the transfer 
equation is provided in Appendix ^ The primed solid angle element dil' is associated with the unprimed direction 
vector n since we wish to reserve dfl for the polar angles of the cluster on the sky. 

We now define a polarization basis to define the Stokes parameters of the radiation incident at the observer from a 
cluster in any direction on the sky. We denote the polarization basis vectors as ei, 62 and leave these unspecified for 
the moment. The Stokes parameters measured in the 61,62 basis at our position due to scattering in the cluster at 
comoving position x are then: 

I{x) + Q{x) = 2Iij{x) ei^eij , 
I{x)-Q{x) = 2Iij{x) e2,ie2,j , 

U{x) = 2hjix) ei,,e2,j . (37) 



Note that we may ignore the Stokes V parameter - it remains zero since no circular polarization is generated by 
Thomson scattering. On substitution of Eqn. H35|l we find 
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Irn 

= -^TcMBiT)Y,aimix) J dn' YUn){ei ■ h){e2 ■ n) . (38) 

Im 

To perform the angular integral we need to expand the integrands in equations (|38|) in spherical harmonics by 
expressing h in polar coordinates. In polar coordinates about the z axis, the Cartesian components of the direction 
vectors are taken to be: 

n ~ {hx,ny,nz) ~ (sin 6' cos (f)' , sin 0' sin (j)' , cos 6') , 

X — (xj:, Xy, i^) = (sin ^^ cos 0, sin sine/), cos 0) . (39) 

We use the following method. Any spherical harmonic can be expanded in terms of the complex quantities (zi, Z2, z^) = 
(sin O'e^"^ , sin 9'e~'^'^ , cos 0') (see for example the discussion of spherical harmonics in ^3). In terms of these functions 
we may write n = ((zi + Z2)/2,i(z2 — zi)/2,Z3). Note that Z2 = Zi, Zg = Z3, Z1Z2 = 1 — z|. It is also convenient to 
define 



Then for instance we have 



e = - iey . (40) 



Zl Z2 

Y ei • e_ + y ei • e+ + Z3 ei • (41) 



The I = 2 spherical harmonics may be written as functions quadratic in the z's as follows: 

Y - J ^(K^ ^ 
-12. n — \ -. — — Zq — — 



'''' - yA^yr^ 2 

VY2^1 = -W— Z1Z3, Y2-l = \ —Z2Z3 

V OTT V OTT 

Y22 = -\f^zl Y2-2^-\f^zl. (42) 
' 4 V 27r ^' ' 4 V 27r 2 ^ ' 

Then we can decompose the integrands into spherical harmonics by expanding the integrands in equation 138|) in 
the z's using expressions like (|41|l and comparing with the expressions for Y2m above. (Note that in performing this 
calculation, it is necessary to use the relation ziZ2 = 1 — z| to eliminate one of the coefficients (zi,Z2,Z3)). We find 
the following manifestly real result for the integrands of Q and U : 

8^ ' 



{e2 ■ hf - [ei ■ nf = \j <3™(*)^2,m(ri) , 



m=-2 
2 



-2(ei.n)(e2-n) = ^ C/,„(«)r2^,„(n) . (43) 

m=-2 

The coefficients Qm,Um appearing in this expression arc the following functions of the arbitrary polarization basis 
vectors chosen by the observer, which in turn are functions of the cluster direction on the sky (so Qm, Um are written 
as functions of the cluster direction vector x, which will become explicit once a polarization basis is chosen): 

Qo{x) = --^[(ei -6^)2 -(62 -6^)2] , 

Qiix) = [(e_ • ei)(e^, • ei) - (e_ • e2)(e3 • 62)] , 

Q2{&) = ^[(e_.e2)2-(e_.ei)'] , (44) 
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and 



Ui{x) = [(e_ • ei)(e^ • 62) + (e_ • e2)(e2 • ei)] , 
U2{x) = --^(e_ • ei)(e_ • 62) . 
Also, we define the quantities with negative m by the relations 

g.,„(*) ^ (-i)"q:„(&) , 
u-„.{x) EE (-i)"c/;;(*) . 

We specialize now to a particular choice of of polarization basis vectors. A suitable choice is 

Bz - l^X 



(45) 



ei 



62 



(46) 



(47) 



where fi ~ x ■ e, ~ cos0, so that the Stokes parameters Q, U are defined with respect to the plane containing the 
axis and the photon direction —x (see Fig.l^J. In this polarization basis the coefficients Qm, Um are 



Qa{x 
Qi{x 
Q2{x 

Ui{x 

U2{X 



— sin^ 
1 „ , 

1 



cost/smWe 



2V3 
, 
i 

i 



1 + cos^ e)e-'^"t 



^/3 



cos f 



(48) 



where 6 is the polar angle between x and the axis and is the azimuthal angle between the projection of x on the 
{ex^ey) plane and the e.^ axis. 

The Stokes parameters may now finally be expressed as a linear combination of the a2m{x), with polar axis e^. 
Angular integration picks out the five coefficients a2m of the primary anisotropy: 

2 

Q{x) = TcPo ^ Qm{x)a2m{x) , 
m=-2 
2 

U{x) = rcPo J2 Um{x)a2ra{x) , (49) 



-—^=Tcmb(.t). Note that this depends on conformal time, but the fractional distortion in the Stokes 



where Pq 

parameters is redshift independent. Also, it turns out that 

4 TStt 



3 V 5 



where 2y/m is the spin- weighted spherical harmonic of spin 2 (see e.g. [T^'l. Thus we may write 

6 [2 ^ 
Q{x)±iU{x) = - — J -tcTcmb{t) ^ ±2>2m(*) a2m(a;) ■ 



(50) 



(51) 



Tn=-2 
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Finally in this section, we quote the following properties of Qm,Um, which are needed in (these are derived 
using the explicit forms in Eqn. (|48|l ): 

2 2 



m — — z m — — z 






2 

E Qm.iXt)U*^{x,) = , 
m=-2 






2 

E (-i)"Q™(*.0O™(*0 

m=-2 


2 

^ 3 


1 

.1-2 


2 

E (-i)™c/™(&,:)f/;(*o 

Tn=-2 


2 

= - cos 26*., 


2 

m=-2 


= . 





(52) 



STATISTICS OF THE CLUSTER POLARIZATION SIGNAL 



Now we consider the information obtainable from measurements of the CMB polarization signal (due to scattering 
of the CMB quadrupole) from galaxy clusters at various redshifts and lines of sight. Assuming that the redshifts 
of each cluster can be obtained, this allows mapping, in principle, of a particular linear combination of aim over 
a significant portion of our past light cone. Galaxy clusters at similar redshifts and on lines of sight separated by 
small angles will produce polarization signals with Stokes parameters which are strongly correlated. Widely separated 
clusters produce uncorrelated signals — and it is the combination of these uncorrelated signals that gets around the 
cosmic variance bound. 

Using Eqn. (|49|l . the two-point correlation function {Q{x)Q{x')) , of the Stokes parameters, as defined in the basis 
Eqn. H47|) . due to two clusters at general comoving positions x,x' with Thomson optical depths tc,Tq is given by 



{Q{x)Q{x')) = tct'cP^ 



2 

E 

m.rn' — — 2 



{x,x') 



(53) 



and similarly for {U{x)U{x')) and {Q{x)U{x')). 

The two-point correlation function {Q{x)Q{x')), for points lying on the same line of sight, is shown in Fig. 01 
which was computed with cosmological parameters Vim = 0.35, f^A — 0.65, n = 1. The solid black curve is the total 
correlation, and the other curves show the contributions from the SW and ISW terms, and their interference term. 
Note that the interference term is negative - its magnitude is shown here. Only the SW contribution has a genuine 
zero crossing. The vertical dashed line indicates the time of (instantaneous) recombination. Note that the SW part 
of the correlation passes through zero at redshift w 10, since at redshifts higher than this x' is in a region of the 
universe separated from the origin by a comoving distance greater than c/{2Hq), and thus the I = 2 correlations die 
off rapidly. 

As a; ^ x', {a2m{x)al^,{x')) C2{T)5jnvi', therefore 

(Qixf) = {U{xf) = \rlP^ C^ir) , 
{Q{x)U{x)) = 0. (54) 
Thus the ensemble average polarization magnitude due to the scattering of the CMB quadrupole is: 



{nixf 



{Q{xf) + {U{xf) 4 2 / Po 



Tcmb(t)^ 



\ TcMB (t ) 



C2{t) 



(55) 



The quadrupole Q2 is conventionally defined by C2 ~ {At: /^)Q'^. Thus the root mean square polarization magnitude 



is given by (recall Pq 



4V107r 



T'cmb(t)) 



(P)^y(n(^ 



V6 
10 



TcQ2 , 



(56) 
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FIG. 4: Normalized magnitude of the two-point correlation function of the Stokes Q parameter, ^P^^'^Tq^\{Q[x)Q{x' ))\/C2{to), 
where x is taken to be at redshift z = Q and x' is a point at redshift z in the plane orthogonal to i. The growth function 
was computed with cosmological parameters = 0.35, = 0.65, n = 1. (Note that the dip in the interference term, which 
occurs at z ~ 3, is not a zero crossing). 



as obtained by In a ACDM model, (P) w 5tc /^K (at zero redshift). so the magnitude of this signal is comparable 
to that of the other SZ polarization effects. 

Now the relations Eqn. (|54|1 suggest estimators C^(t) and C;^(t) of C2(t), given the measured Stokes parameters 
of N clusters at the same redshift z{t) on lines of sight Xi with optical depths Tc.i [i ~ I,- - -jN), which beat the 
cosmic variance limit: 

i 

d^r) ^ lp,-'J2^.U{x,r . (57) 

i 

Note that these are non-optimal, coordinate-dependent estimators. In future work it would be interesting to construct 
optimal estimators of the power at a given scale, given the cluster polarization data, but we do not attempt that here. 
For the mean of these estimators to equal C2(t), the weights Wi must be chosen to satisfy 

Y.W,T^, = 1. (58) 

i 

We will consider the simple choice 

= ■ (59) 



14 




Note that n = would be a bad choice since it gives more weight to chisters producing weaker signals, reducing the 
signal to noise. A better choice is the uniform weighting n = 2. 

The cosmic variance limit on these estimators is determined by the variances (with X, X' indicating cither Q or U) 



Varaf(r) = | V ^ W^. M^, (X(a;0'X(a;,) 



(60) 



The sum over i,j may be broken into a contribution from clusters at the same location, Vari(X), and a contribution 
from clusters at separate locations, 'Va.r2{X,X): 



Var C2 (r) = Vari(X) + Var2(X, X) , 



(61) 



where 



1 + (1- -sin220,)' 



i 

Vari(C/) = Po-^C2iTf ^ W^r^^, [l + cos^ 29,] , 



(62) 



and 



Var2(X,X') = - W^Wj^h^' 



2 



h 3>i 



^ ^ ^yj^l ("^2 )^m2 ) ^2mi 2m2 ("^^ ; '^J ) 



+ 



^ {-l)'^'^X^^{Xi)Xll^{x.j) C2,m2,m2{x„Xj) 



(63) 



Here we used the 4-point correlation function from Eqn. (jSJ, the relations in equations H52|) . and the fact that 
C2m2m'{xi,Xi) = (5mm' C2 (t) ■ Notc that the variances are functions of the angular positions of the clusters on the 
sky, which is due to the specific choice of polarization basis for the Stokes parameters. 

If all of the correlations C2m2.m' {xj , Xj ) with i ^ j vanish, then only the first term on the right hand side of equation 
()61|l remains. If however the cluster positions Xi are close enough that the these correlations approach C2(t), then 
the O(N^) terms in the YslT2{X, X') terms combine to swamp the first terms. Thus in order to beat cosmic variance 
by a factor of 0(7V-i/2) we need sets of clusters which are mutually uncorrclatcd (as pointed out in a qualitative 
discussion by P|). The number of uncorrclatcd regions available increases as the redshift increases, since the comoving 
region surrounding each cluster outside of which the polarization is approximately uncorrelated with that produced 
by the cluster is smaller at higher redshift. This is because smaller comoving scales contribute to the I = 2 harmonic 
of the CMB on the sky at higher redshift, and similarly the same comoving scale maps into different angular scales, 
as illustrated in Fig. [S] C2(t < tq) depends on fluctuations of smaller scale than C2(to). 

At very low redshift, any cluster will be correlated with any other, and we get back the usual cosmic variance 
constraint. In other words, we can beat the cosmic variance on C2(t < tq), but not on C2(to), today's quadrupole. 
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To demonstrate the reduction in cosmic variance, we first combine the Q and U measurements to obtain an improved 
estimator of C2(t). Taking a Hncar combination of C^, yields an improved estimator, 

(r) ^ aC^ir) + (1 - a)d'^ [t) , (64) 

with Q < ot < 1. This has variance 

\&tC^{t) = a^N&T [t) + {\ - afN&T {t) 

+2a(l - a) Cov{d§,d^) . (65) 

It is easy to see that the covariance Cov(C^,C2^) is zero because of the relations in equations H52(l . which imply 
Var2((5, U) = 0. Thus the optimum value of a is trivially a = 1/2. 

With these expressions we can compute the variance of our estimator for C2 given a set of cluster positions and 
optical depths {jCj}, {tci}- For simplicity, we will compute here only the variances for sets of clusters which lie 
on a given redshift shell, distributed in random directions. In the left panel of Fig. we show the variance of the 
estimator for C2 obtained from hypothetical cluster polarization measurements from sets of 10 and 100 clusters, all 
on the same shell of redshift z, distributed in random directions on the sky. (We assume the cosmological parameters 
= 0.35, JIa = 0.65, power spectral index n = 1). In the left panel, we show the variance of the estimator for C2, 
as a fraction of C2 at z = 0, from separate sets of clusters confined to redshift shells as a function of the redshift of 
the shell. 

For simplicity we have assumed that all clusters have the same optical depth tq, so the weights given in equation 
(|59|l are all equal, Wi = t^^'^ /N. The filled squares indicate the square root of Var C2 , the variance of the average 
estimator for C2, (C^ +C^)/2, for 10 and 100 clusters as indicated. This is expressed as a fraction of the Sachs- Wolfe 
contribution to C2, which is independent of redshift for n = 1. The horizontal dashed line at 1/2/5 ~ 0.63 indicates 
the cosmic variance limit on a single CMB sky given in Eqn. @. (Note that the ISW contribution is included in 
these calculations, but omitting it leads to a difference of less than a few percent in the curves in Fig. |SJ). As z — > 0, 
the estimator variance slightly exceeds the cosmic variance limit given a single CMB sky — clusters with overlapping 
correlation spheres as z — s- are no more useful than a direct measurement of the quadrupole on our sky by e.g. 
WMAP. In fact our estimator of C2 is worse than a direct measure — but an optimal estimator could be constructed 
which would yield all of the a2m at z « from cluster measurements at various points on the sky. As z — * 00, 
the estimator variance goes as ^ 1/VN (as z 00, the variance approaches 0.923/ Vn. This may be derived by 
averaging the estimator variances in Eqn. ()62(l over angles in the N —> 00 limit). 

Thus if we have enough clusters, we can measure C2(z) at increasingly high accuracy as the redshift shell increases. 
However the fact that C2(z > 0) can be measured more accurately than C2(z = 0) does not necessarily mean that 
we can reconstruct the primordial power more accurately than we could using only say the WMAP multipolcs, unless 
the variance in our estimator is less than the variance in the CMB multipole which probes the power at the scale 
corresponding to the quadrupole at redshift z. 

In the right panel of Fig. Elwe show the variance of our estimator of C2, except divided through by the cosmic 
variance in the CMB multipole Ci^,.f , where leB is the Z-scale on our sky corresponding to the I = 2 scale on the redshift 
shell at z > 0. To be explicit, we are plotting 

\/varC72 

, ^ ■ (66) 

v/2/(2/eff + 1) 

This is the key plot, since if this ratio is less than unity, we have shown that this method can improve on the cosmic 
variance limits inherent in the existing CMB data, as discussed in the introduction. A reasonable approximation for 
Ics is to take the real I value corresponding to the comoving radius of the last scattering sphere at the redshift z, 
namely Zoff = 2/(1 — x/(''"o — ^r)) where x(z) is the comoving distance of the shell. Note that there is no one length 
scale or (z = 0) harmonic which corresponds to the (z > 0) quadrupole, since a single angular harmonic mode has 
power spread over a range of scales. In Appcndix^we check the approximation for l^s given above by considering the 
fc-space window function of the quadrupole, and find that a slight modification to the formula for /cff given above gives 
a better approximation. We use this modified formula in our subsequent calculations. We opt to compute l^ff as a real 
number and interpolate between the Cj's at the integer values of I bracketing Zoff, to obtain the corresponding cosmic 
variance. This will give a reasonably good approximation to the accuracy achievable in power spectrum estimation. 

A rough estimate of the maximum number A'niax of clusters with uncorrelated signals available on a given redshift 
shell is given by dividing the area of the shell by the area of a circle with radius equal to that of the last scattering 
sphere at that redshift: 



iVmax = 1 + 4xV(7-o - Tr - xf , 



(67) 
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FIG. 6: Reduction in cosmic variance with clusters distributed on various redshift shells in random directions. 



where is the comoving distance to the shell. The first term is added to allow for the fact that at least one 
independent cluster signal is available at low redshift. An approximation to the minimum estimator variance achievable 
on a given redshift shell is then given by: 

V275 C.(.^_0)_ ^g^^ 



An approximate lower bound on the ratio of the square root of the estimator variance and the cosmic variance in Ci^j^ 
is thus given by: 



"-^/cft V -'vmax 

This is shown by the solid curve in the right panel of Fig. Since this curve exceeds 1 for all values of z, this suggests 
that the cosmic variance is not reducible below the WMAP levels with the cluster technique even without computing 
the detailed correlations. 

In the right panel of Fig. El the black squares joined by dashed lines show the ratio (|66|l computed with the full 
expression for the estimator variance, plotted for sets of 5, 25 and 100 clusters as indicated. Clearly there is no 
reduction in variance below the WMAP cosmic variance at any redshift with up to 100 clusters, as expected from 
our rough analysis. The three curves are barely distinguishable up to z « 10, reflecting the fact that at low redshift 
increasing the number of clusters does not enhance the reduction in cosmic variance since the polarization signals are 
correlated. Increasing the number of clusters further will not lead to significant improvement, at least with signals 
from z < 10. For z > 10, the = 25 and N = 100 cases begin to show a drop in the ratio, as the redshift shell 
expands (and last scattering spheres around each cluster shrink) to the point at which the signals from each cluster 
are uncorrelated. As z increases eventually the ratio in the A^ = 25 and A^ = 100 cases begins to grow again, since 
Icff is increasing as the last scattering spheres shrink, and so the cosmic variance in Cj^jj. is decreasing. 

Unfortunately it seems that the cluster technique, at least using the estimator we propose, is not competitive with 
the CMB data at z = for estimation of the power at a given scale. In future work it would be interesting to attempt 
to construct a better estimator — for example combining signals from sets of clusters distributed at arbitrary redshift 
rather than on redshift shells would presumably lead to some improvement. However we emphasize that this technique 
is still a useful probe of the time evolution of C2, and possibly the ISW effect. 
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VI. DISCUSSION 

We have developed a statistical theory of the part of the polarization signal in the CMB in the direction of galaxy 
clusters produced by scattering of the CMB temperature quadrupole. We have shown explicitly that it is possible 
to use the indirect information about the last scattering surfaces of distant observers contained in these polarization 
measurements to constrain the / = 2 angular power spectrum harmonic of the CMB, C2, as a function of redshift with 
greater statistical accuracy. We also showed however that it not possible to use the cluster polarization measurements 
to probe the power on a given scale with higher accuracy than the limits imposed by cosmic variance on the single 
sky CMB data. Thus power spectrum estimation cannot be improved using this method. 

But we believe the cluster method is still of considerable value though, since it serves as a probe of the physical 
mechanism which might have suppressed the quadrupole. It has been noted that the quadrupole seems to be lower 
than might be expected due to a statistical fluctuation alone 0|. The quadrupole may be anomalously low, that 
is lower than the standard models predict, for various reasons: a cutoff at large scales in the fluctuation power 
spectrum, different-than-cxpected behavior of the transfer function at large scales, or effects associated with the 
large-scale topology or geometry of the universe. There could conceivably be other explanations, but these are the 
simplest. Signatures of these effects could be present in the time evolution of C2. For example if the quadrupole is 
low because of some topological suppression, we should see a rise in C2 with redshift as the scale probed falls below 
the local quadrupole scale (but if the standard models are correct, there would be no such rise, at least in the n = 1 
case). 

One might be worried that since the accuracy of the cluster measurement of C'2 is not higher than the accuracy of 
the measurement of the corresponding WMAP harmonic, as shown is Jvl the cluster method may not provide any 
more information that that already contained in the WMAP data. However the cluster technique yields information 
which is not obtainable with WMAP (about perturbations on last scattering surfaces different to our own) and is thus 
complementary. In future work it would be desirable to perform a full analysis of the feasibility of using the cluster 
method to detect the effects of topology (or other quadrupole suppression mechanisms) in the time evolution of C2 , 
and a comparison with what we can already learn from the WMAP data. 

We also showed that in the standard ACDM model the ISW effect produces a small (~ 2%) bump in CMB harmonic 
C2(t), which is swamped by the high cosmic variance at low redshift even using large numbers of clusters. However 
the ISW effect produces a significant feature in the two-point correlation function of the Stokes parameters which 
might be detectable. Detection of the ISW effect would provide additional information about the acceleration of 
the universe and the dark energy. We also note that this method is a rather sensitive probe of deviations from the 
scale- invariant power spectrum, because if n 7^ 1 then the Sachs- Wolfe contribution to C2 (t) either grows or decays 
rapidly with conformal time. 

The procedure we have outlined is something of an idealization. We assumed that the polarization signal induced 
by the quadrupole is obtainable from many clusters at the same redshift, and we ignored noise and contamination of 
the signal. In practice there are contaminating polarization signals from the kinetic and thermal Sunyaev-Zeldovich 
effects and distortions in the polarization field due to lensing |2^. Clearly separating the quadrupole signal 
from the other contaminants would be a major experimental challenge. However the signal to noise may be increased 
to some extent by combining signals from clusters located at similar directions and redshifts, since the signal from 
sufficiently nearby clusters is strongly correlated. It remains to be seen if this technique will be a practically useful 
cosmological probe. 
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APPENDIX A: TRANSFER EQUATION FOR POLARIZATION MATRIX 

We found it convenient to write the transfer equation for generation of polarization by Thomson scattering in matrix 
form. This approach is similar to the "density matrix" formalism for polarized radiation transfer 0, . 

The transfer equation is usually written in terms of the Stokes parameters, which are time averages of quadratic 
products of the complex amplitudes of the electric field components Ei of the electromagnetic wave. The polarization 
state and intensity of a beam of light propagating in the z-direction is characterized completely by the 2x2 Hermitian 
matrix (EiE*), with (z, j) G {x,y} (the brackets denote a time average). An obvious generalization is to allow 
to become Cartesian tensor indices and to run over all of {x,y, z}. We obtain a 3 x 3 polarization matrix associated 
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with photon direction n: 

Q,j{n)^{E,E*) . {i,j)€{x,y,z} (Al) 

With this matrix wc arc no longer restricted to considering a beam propagating along a coordinate axis. For a given 
photon direction n, the electric fields are transverse to n, implying 

n'gy(n) = 0. (A2) 

Consider now a superposition of beams of various directions n and momenta p ~ pn, where p is the frequency (we 
set c = h = 1 for convenience). In this case, Qij cannot be considered a function of the photon momenta, but an 
intensity matrix can be defined associated with each photon direction and frequency. Recall from the definition of 
specific intensity that the energy density is given by = J I{p) dpdil, where I{p) is the specific intensity and dfl is 
the solid angle element associated with the photon direction n. Similarly we can express each element of the matrix 
Qij in terms of specific intensity matrices : 

Qtj = J dpdVthj{p) . (A3) 

Now the Stokes parameters are defined with respect to a particular choice of "polarization basis" . This is a pair of 
mutually orthogonal unit vectors e^^\e^^\ both orthogonal to the beam direction. The Stokes parameters are given 
in terms of lij and the polarization basis vectors as: 

'^j -'-V 1 
^(2)^(2) J.. 

,(l)g(2)j.. 
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+ Q 
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-Q 
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>riV 
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u- 


-iV 



^^I. ■ (A4) 
In the case of a beam propagating in the z-direction for example, we have, choosing polarization basis vectors 



I + Q U + iV 
~iV I-Q 




h,^T;\ U ~iV I^Q Q \ . (A5) 



The matrix is non-zero only in the two dimensional subspace spanned by e^^-*, e^^-'. 

We need to construct the matrix of an unpolarized beam propagating in a general direction n. The only quantities 
available to form this matrix are the intensity /, the components of the direction vector n, and the Kronecker delta 
5ij. The matrix must therefore be of the form: 

lij (n) = ASij + BniUj . (A6) 

The matrix of an unpolarized beam propagating in the z-direction is obviously 

1 
1 




/y = ^ I 1 I . (A7) 



Comparing this with the form of Eqn. (|A6p for the special case rii = 5iz, we see that A = —B = J/2. Thus the matrix 
of an unpolarized beam in a general direction n is 

= ^((5y - n^Uj) . (AS) 
The polarization magnitude of the beam described by a general matrix lij is given by 

= 2Tr[j2]/Tr[/] - 1 , (A9) 
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which may be readily checked with the matrix (|A5|) . 

We now derive the equation for the time evolution of the polarization matrix due to Thomson scattering from a 
distribution of stationary electrons. For a completely linearly polarized beam, Q^^£*£j is the time-average energy 
density for electromagnetic radiation of polarization e, where e • e = 1. Consider a completely polarized beam with 
polarization vector e and momentum p = pn incident upon an electron at rest The polarization matrix of the incident 
beam is Qij = Qeie* where cQ is the incident flux (we choose units such that c = 1). Normalization of the polarization 
vector implies Q = Qije*ej. In the Thomson limit, in which the electron recoil is negligible, the differential cross 
section for Thomson scattering of a beam into final momentum p' = p'n' and polarization e' is [2l] | 

where dfl' is the solid angle element associated with the scattered photon direction n' . Thus the power per unit solid 
angle in the scattered beam is 

We may also write Q\e* • e'p ~ Qije'*t'y 

Next consider a gas of electrons all at rest with number density rig. We can ignore the thermal motion of the 
electrons here since the correction to the polarization due to a finite electron temperature Tg will be down by a factor 
of 6e = kBTe/m^c^ . At cluster temperatures 9,, « 0.01, so we may assume stationary electrons. Assuming incoherent 
scattering, multiplying Eqn. (|A11|) by riedfl' converts scattered power per electron to the rate of change of energy 
density in final polarization state e': 



dQ'ij 1^ I _ 3(Tt 
IT ~ "8^ 



e'*e', = ^n,Q,je'*e'j dn' , (A12) 



where Q'^^ is the polarization matrix of the scattered beam. Using equation (|A3|I , setting p ^ p' since we are working 
in the Thomson limit, and assuming the incident beam is monochromatic, we find 



- "s^^^T-'e H J dn Uj , (A13) 

where d^ is the solid angle element associated with the incident photon direction n. Note that this equation is valid 
now for an arbitrary incident radiation field. 

We cannot now just remove the polarization factors and conclude dl[j oc / hj because the polarization of the 
incoming wave does not lie in the same plane as the polarization of the scattered wave. For a given outgoing momentum 
p', the outgoing polarization is a linear combination of two basis vectors e'^ and e'2 (orthonormal and orthogonal to 
the photon momentum p'). Thus, lijt'* e'j projects out of the incoming matrix lij only those components lying in the 
e'j^-ej plane. This projection is equivalent to first projecting Jjj with e\ (g) e\ + ® e'j. This is equivalent to projecting 
out the unphysical components by acting with the matrix 

P,j(n') ■ (A14) 

Projecting the final polarization vector with Pij{n') does not change it: Pij{n')e'j = ej. It follows that Iij{n)e'*ej = 
Iki{n)Pik{n')Pji{n')e'*ej. Now it is safe to remove the outgoing polarization vectors from Eqn. (jAlSp . 
We conclude that, for any initial and final polarizations, 

^Il2pl ^ ^n.P^kin') P,i{n') f dil hi{n) . (A15) 

at OTT J 

This is the equation we need. The projection tensors are easy to understand: the scattered matrix is simply propor- 
tional to the incident matrix after the unphysical polarization components (those proportional to n') are eliminated. 

If the integration time is sufficiently short, we may replace rieCT dt with the optical depth to Thomson scattering, 
re- Then we have for the polarization matrix of the scattered radiation field: 

4, = ^P^k (n') Pji (n') j dn hi (n) . (A16) 

If the incoming radiation field is assumed to be unpolarized, then the incident matrix has the form of Eqn. I)A8|1 and 
we may set /fe;(n) = lin) Pki{n)/2. Thus 

/^ (n') = ^P^k{n') P,i{n') j dn I{n) Pui{n) . {All) 
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FIG. 7; The lower, middle, and upper sets of points indicate the wavenumbers at the 20th, 50th and 80th percentile respectively 
of the window function for the I — 2 harmonic, as a function of redshift. This gives the range of wavenumbers contributing to 
the CMB quadrupole at a given redshift. The lower line shows the estimated wavenumber fcofi corresponding to the radius of 
the last scattering sphere. The upper line shows k'^g = 3.12 fcoft , which better matches the median of the window function. 

This is the form of the transfer equation used in Eqn. (|36|l . 

Note that this is not the full transfer equation, since we have ignored the effect of scattering out of the beam 
direction. But here we are only interested in the polarization generated by scattering into the observation direction, 
and the loss of photons from the beam only affects this at 0{tq). 

Assuming a blackbody spectrum of incident photons, this transfer equation also holds for the brightness temperature 
polarization matrix (since Thomson scattering docs not change the photon frequency, and there is no Doppler shift 
because we have assumed the electrons are stationary, thus the scattered radiation field also has a blackbody spectrum). 
We may suppress the frequency dependence of the (brightness temperature) polarization matrix and associated Stokes 
parameters of the scattered photon also, since there is no energy transfer. 



The I = 2 harmonic of the CMB anisotropy at redshift z = contains contributions from a broad range of 
wavenumbers of order c/Hq. At higher redshifts the CMB quadrupole probes the potential on a smaller last scattering 
surface than at redshift z ~ 0. In Wl we need to associate the quadrupole at z > to an effective angular harmonic 
mode /eff probing roughly the same angular scales on the z = sky. This is most easily done by finding the comoving 
wavenumber fceff which best approximates the scale probed by the quadrupole on the z > last scattering sphere, and 
setting Icff = 2 fces/fco where fcg is an estimate of the scale probe by the quadrupole at z = 0. We take /cq = l/(''o — tv) 
(we have set c = 1). 

An approximation to the effective length scale probed by the CMB quadrupole at a given redshift is given by the 
radius of the last scattering sphere at that redshift, r = (tq — r^) — x (where x(z) is the comoving distance to redshift 
z). With fcoff = i/r this gives ks = 2/(1 - x/ (^o - ^r))- 

To check this approximation, there is a useful method which yields an estimate of the range of wavenumbers k 
which contribute to each / mode. The procedure is to first find the following window function: 



APPENDIX B: EFFECTIVE SCALE OF THE z>0 QUADRUPOLE 



P{k\T,l) 



(47r) 



Ci{r) 




(Bl) 



where the transfer function is computed with some assumed model. In the I = 2 case the Sachs- Wolfe transfer 
function is sufficient until we consider rather high redshifts. The k scale probed by this harmonic is then taken to be 
the median ki of the window function. The range of wavenumbers contributing to the I harmonic at a given redshift 
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is indicated by the 20th to the 80th percentile of the window function (which corresponds to full-width-half-maximum 
in the Gaussian case). The I = 2 case is shown by the points in Fig. \7\ For comparison the estimate kgg is shown by 
the lower line. In fact we obtain a much better match to the median window function estimate with the wavenumber 
k'^ff = 3.12 kcs, as shown by the upper line. We therefore use the following modified approximation for the effective 
angular harmonic mode in the calculations in ^ I'^ff = 6.24/(1 — ^/(to — Tr))- 
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